
#include "stdafx.h"
#include <stdexcept>
#include <iostream>
#include <map>
#include <ctime>
#include <cstdlib>
#include "string.h"
#include "blinalg.hpp"
#include "superelement.hpp"
#include "mesh.hpp"
#include "gauss.hpp"
#include "Intergration.hpp"
#include "functions.hpp"

FullMatrix<double> computeStiffnessMatrix(vector<FElement> ip_elems);
EVector<double> computeRightHandSide(vector<FElement> ip_elems);
double g(double ip_x);
//double app_solution(EVector<double> ip_c, vector<FElement> ip_elems, double ip_x);
EVector<double> giaiHePT(FullMatrix<double> B, EVector<double> F);
double computeL2Error(vector<FElement> ip_elems,EVector<double> ip_c);
double computeL2Error(vector<FElement> ip_elems, EVector<double> R, vector<FElement> m_elems);
double computeEError(vector<FElement> ip_elems,EVector<double> ip_c);
